Comments on the Cooperstock-Tieu Galaxy Model 
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ABSTRACT 



ON ! The recently proposed Cooperstock-Tieu galaxy model claims to explain the 

flat rotation curves without dark matter. The purpose of this note is to show 
that this model is internally inconsistent and thus cannot be considered a valid 
Q\ \ solution. Moreover, by making the solution consistent the ability to explain the 

flat rotation curves is lost. 
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c/2 . 1. Introduction 

Cooperstock and Tieu model a galaxy in general relativity as an axially symmetric, 
pressure free dust cloud with metric 

ds 2 = -e w (cdt - Ndcp) 2 + r 2 e~ w d<p 2 + e v ~ w (dr 2 + dz 2 ) (1) 

where w, v, and iV are functions only of r and z, thus the metric is stationary (Cooperstock 
& Tieu 2005). They further assume their coordinates to be co-moving with the galactic dust, 
thus 

v? = e - w/2 5? (2) 

is the four- velocity. 

By performing a local diagonalization of the metric they obtain a relation between the 
metric and local angular velocity as 

cNe w cN 
= ~ (3) 
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where the approximate value is appropriate to the order they consider. The tt and tip Einstein 
equations are respectively 

r-^N^Nl) = ^ (4a) 
N rr -^N r + N zz = 0. (4b) 
The second equation is equivalent to the formula 

V 2 $ = (5) 

where they have defined 1 




Thus the observed rotation curve becomes a boundary condition for the solution to 
Laplace's equation (5) which they take in the form 

<b = Y,C n e- k ^Mk n r) (7) 

n 

where the k n are chosen for orthogonality over the radius of the galaxy Once N is found 
by fitting this function to the obsersved rotation curve, they derive the density by (4a) and 
in this way they obtain an excellent fit to the data while obtaining a density profile that 
accords with observation. 

However, it has been pointed out (Korzynski 2005; Vogt & Letelier 2005) that, since the 
solution depends on \z\, equation (4b) is not satisfied, but rather yields a singular contribu- 
tion to the z = plane, which has the properties of an exotic form of matter. It may be 
wondered whether this singular disk can be removed by choosing a different solution form or 
by increasing the complexity of the model. However, in the following analysis we will show 
that this is not possible and that model fails to accord with general relativity. 



2. Analysis 

Assuming this form of the metric, and without making any approximations, the scalar 
of volume expansion = u 1 * vanishes (a semicolon denotes covariant differentiation and a 
comma denotes partial differentiation). Defining the space-projection tensor 

/V = 9v» + u v u » ( 8 ) 



It should be noted that this $ is not the Newtonian gravitational potential. 
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the shear tensor is given by 2 

<V = u {c;P)h a li h p v (9) 
and this vanishes as well 3 . However, the vorticity tensor, given by 

c<V = u [a ^]h a ^ u (10) 

does not vanish, but has nonzero components 

u w = -u rv = ±e w/2 N r (11a) 
u<r» = -u*p = \e w ' 2 N tZ . (lib) 

Though the matter in this model does indeed rotate, the rotation is rigid and thus cannot 
characterize a galaxy which is differentially rotating. It should be emphasized that since 
is a tensor, this cannot be a coordinate effect. 

With the present solution Raychaudhuri's equation (Ciufolini & Wheeler 1955) simplifies 

to 

= R tt /gu (12) 

and in fact reduces to the condition V 2 w = 0. This condition is demanded in Cooperstock 
& Tieu (2005) on the grounds that the geodesic equation be satisfied. This is equivalent to 
saying that the geodesies must be circular orbits about the z-axis, which should not hold in 
general. Orbits in the z = plane should indeed be azimuthal, but we cannot expect this 
behavior off that plane. That Raychaudhuri's equation demands this condition again reveals 
the rigidity of the rotation 4 . 

Now, if we seek solutions to (4b) that are symmetric about the plane, and singularity 
free, then must require iV to be independent 5 of z. Thus (4b) has the trivial solutions 

N = A or N = Br 2 (13) 

where A and B are constants. The first solution leads to zero density and the second 
to a constant density under rigid rotation, according to equations (3) and (4a). Thus it 



2 {pi,i') and \jjb,v] denote symmetrization and antisymmetrizaion with respect to the enclosed indices, 
respectively. 

3 The vanishing of shear was also found by Bonner in his solution (Bonner 1977). 

4 Actually, any co- moving coordinate system requires g tt — —1 as the coordinate points are in free fall 
and thus keep proper time (Weinberg 1972), which here requires w = 0. 

5 We could choose cosh(z), but this would lead to an exponentially increasing matter density. 
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appears that the physical origin of the singularity is in attempting to describe, in co-moving 
coordinates, a non-rigidly rotating dust cloud, which the metric (1) cannot. 

Next, in order to solve the Einstein Equations Cooperstock and Tieu perform an ex- 
pansion of the metric in \/G and conclude that the functions w and v are of second order, 
but the function N, which couples to the rotation, is of first order. Strictly speaking, this 
expansion is not well defined as the expansion parameter has dimensions. We can form the 
dimensionless parameter 

A = V^ (14) 

where M is some characteristic mass of the system and L some characteristic length (for 
example, the mass and radius of the galactic core). We then compare equations (4b) and 
(3): 

AV(A;.A:) = ^ Gp 



c 2 



A— 

r 



where the order has been shown explicitly. Substituting v for N in the first equation yields 
the relation 

v = O (V^GVlL 2 ) (15) 
where we have taken derivatives to be of order 1/L. This can be rewritten as 



v = (16) 



M 

and since the quantity under the square root is of order unity, we have 



which is expected from Newtonian theory and is the basis of the PPN expansion. 

Now, suppose we choose a coordinate system in which the galactic dust has coordinate 
velocity cu/c, so that the stress-energy tensor has the form 



/ 1 - \ 

c 



oc 



c \cJ 



(18) 
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in the t9?-subspace. The t<^-Einstein equation then has the form 

G t(p 



c 2 c 



(19) 



up to a factor of order unity due to the constraint condition -u M w M = —1. Thus, the right 
hand side of this equation begins at order \ 3 /L 2 according to (17), whereas the left hand side 
is of order \/L 2 according to (4b), since N is assumed to be of order A. Thus the assumption 
that N is of first order is inconsistent, while consistency requires that N be of at least third 
order 6 . 

Moreover, given the form of the stress-energy tensor above, suppose we make a global 
transformation <p — > <p + u>(r, z)t to the co-moving frame, so that all components of the new 
stress-energy tensor vanish except the density. The new metric will have the same form as 
the old metric, but for the differential dtp' we have 



dip' = ujdt + dip + t{uj jT dr + w tZ dz) 



(20) 



which necessarily introduces time-dependence into the new metric unless uj is spatially- 
independent, that is, unless the rotation of the matter is rigid. Thus, contrary to the 
assumption of Coopertock and Tieu, the metric (1) cannot both be co-moving and time- 
independent. This accords with zero value of the shear tensor above. 

It can be seen in the following figure that a co-moving metric of a differentially rotating 
system is time-dependent and possesses shear. Here, r-coordinate lines "twist" up in time 
relative to observers at spatial infinity. <p — ujt has been plotted, where u(r) is the fit to the 
Milky Way from Cooperstock & Tieu (2005). 




Fig. 1. — The twisting of r-coordinate lines for various ip at two different times as seen by 
observers at spatial infinity. 



6 This will be demonstrated more explicitly in the next section. 
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3. Non-co-moving Expansion 

In this section we carry out an expansion of the metric (1) in a system of reference in 
which the galactic dust has coordinate velocity 

U - = *t = W - (21) 
u° dct c y ' 

which is physically the angular velocity measured by observers at spatial infinity (Bardeen & 
Wagoner 1971). In this direct approach we will show that the angular momentum coupling 
is too weak to account for the flat rotation curves. 

We expand the metric (1) as 

9tt = -l-w-w+0(X 6 ) (22a) 

g tv = N + wN + N +C(A 5 ) (22b) 

g w = r 2 - N 2 -r 2 w +0(\ 4 ) (22c) 

9rr = g ZZ = l+W-l+0{\ A ) (22d) 

where the over-script indicates the order of the term in A, which is the same as defined 
in (14). The presence of only even terms in the and odd in g tip are for the appropriate 
behavior under time-reversal (Weinberg 1972). 

The constraint condition = — 1 requires that 

2 

u l = (g tt + 2^ gtv + ^g w y 1/2 (23) 
and thus the stress energy tensor has the expansion 

/ 2 2 \ 

T tt = pc 2 (l-^+L^-) +0(\ 4 ) (24a) 



T tLp = pc 2 - + 0(X 3 ) (24b) 

c 

L) 2 

T w = pc 2 — + C(A 4 ) (24c) 

As seen in the previous section, 87rG/c 2 increases the order by two so that the right-hand 
side of the Einstein equations are of second order and higher, thus the tip equation to first 
order is i 1 i i 

N,rr -- N, r + N ,zz= (25) 
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i 

so that the lowest order term of N is sourceless. We are free to choose N however we wish 

1 

to make the solution simplest, thus we set N= 0. With this selection the Einstein equations 
through third order become 

-V 2 w +\ (v, zz + u >rr ) = (26a) 
v 7r = (26b) 

V,zz + V,rr = (26c) 

v,z = (26d) 

i 9 ( 3 13 3 \ 8nGpru 
-\r- 2 [N,rr— r N,r + N, zz j = ■ (26e) 

2 

We see that v must be a constant so that we have 

V 2 ^ = (27a) 
3 13 3 lQ7rGpr 3 u , „, . 

N,rr~-N,r + N, zz = . (27b) 

Thus we see that the coupling to the angular momentum is of third order, there is no longer 
a nonlinear term in the mass density 7 , and we can identify 

&= - 2 4 (28) 
c z 

with the Newtonian gravitational potential. 

If we analyze circular orbits on the plane the geodesic equation demands that 

r£ + 2r^ + r^ = o (29) 

which to third order can be written for (j, — r as 

V 2 ^ A V 2 \ VC 3 

- = -$ r 1 + - - - N, r 30 

r \ c 2 J r 

which is recognized as the Newtonian centripetal equation plus second order corrections 8 . 
Thus the matter essentially moves according to the predictions of Newtonian gravitation 



7 Even when G u is written to fourth order the nonlinearity due to N is not present. 
8 The presence of the c in the last term effectively lowers the order by one. 
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with corrections of order v 2 /c 2 , which cannot account for flattening of the rotation curves 
without extra non-luminous matter. 

Finally, we can compute the shear tensor, which has the non-zero components 

(m*)V 

(V) 3 r 2 

otz = o z t = (31b) 



2c 
(w') 3 r 2 



0>r = &rtp = ~ W,r (3lCJ 



&<pz = o ztp = — — — (31d) 



all of which vanish exactly when uj is constant. 



4. Conclusion 

It has been shown that the Cooperstock-Tieu galaxy model is inconsistent as a general 
relativistic model and that a proper model fails to account for the flatness of the rotation 
curves without the dark-matter hypothesis. This failure is due to the weakness of the metric 
coupling to the angular momentum of the galaxy. 

However, the flat rotation curves seem to imply a large inertial induction effect, where 
the rotating inner matter boosts the rotation of the outer matter, leveling off the rotation 
curve, which is what the Cooperstock-Tieu model attempts to describe within general rel- 
ativity. Since their solution predicts a matter density well within visible limits it is quite 
possible that their solution represents an alternative, more Machian, gravitational theory 
where inertial induction effects are much larger than in General Relativity. 

The author is greatful to Prof. R. Gilmore for his encouragement and many helpful 
discussions on this project. 
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